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NATIONAL ADVISORY COMMITTEE FOR AEROMUTICS 

ADVANCE RESTRICTED REPORT 

THE TU0-Ip4ENSI0NAL INCOMPRESSIBLE POTENTIAL FLOW 
OVER CORRUGATED AND DISTORTED INFINITE SURFACES 


By W. Perl and L. J. Green 


SUMMARY 

The two-dimensional incompressible potential flow over corru- 
gations and humps of arbitrary shape is derived by conformal trans- 
formation. The results are compared with those obtained by the 
methods of thin-airfoil theory. Some discussion is included of the 
flow over bumps that protrude both inward and outward from a wall. 


INTRODUCTION 

Analyses of the effects of local surface distortions on the drag 
and the critical speed of airfoils usually begin with a consideration 
of the two-dimensional incompressible potential flow over a surface 
having these distortions. In reference 1, for example , the well- 
known approximate methods of thin-airfoil theory are applied to the 
calculation of the velocity distribution over periodic corrugations 
and isolated bumps of sinusoidal shape. 

In this paper the ideal flow past such shapes is derived by more 
exact conformal mapping methods of ¥. Perl of the NACA staff, particu- 
larly inasmuch as the numerical application of these methods is almost 
as simple as that of the approximate methods of thin-airfoil theory. 
The res’ll t s are compared with those obtained by thin-airfoil theory. 
Some incidental discussion of the conditions at a cusped edge and 
of the mapping of bumps that extend both inward and outward from a 
wall Is also given. 

The analysis in this paper was begun at Langley Memorial Aero- 
nautical Laboratory and completed at the Aircraft Engine Research 
Laboratory of the NACA at Cleveland, Ohio. 
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THE FLOW ora A CORRUGATED SURFACE 

Consider a corrugated surface, all the cross sections of which 
parallel to a fixed plane are the same, having infinite length and 
arbitrary hut periodic 3hape. It is desired to find the velocity 
distribution produced along the surface by an ideal incompressible 
fluid moving parallel to this plane. The free-stream velocity suf- 
ficiently far from the corrugation is assumed to be constant, par- 
allel to the axis of periodicity of the cross section, and of 
magnitude unity. 

The problem is solved by finding the conformal transformation 
between points of the corrugation (actually the cross section), 
taken as periodic about the 6-axis of a z-plane, and points of a 
straight line, taken as the 0-axis of a £ -plane (fig. l). The 
Cartesian mapping function (CMF), which relates conformally corre- 
sponding pairs of points in the two planes, is defined as the vector 
difference z - £ between such pairs of points. 

Thus 


z = + i 0 

£ E Wq + i0 

z - £ £ Q - i€ = (i|r - \J/ 0 ) - i (0 - 0) 
The various quantities are defined in figure 1. 


( 1 ) 


The CMF z - £ can be regarded as a function that is regular 
everywhere outside a circle by virtue of the transformations 


z - log p 1 
£ = log p 

in which the coordinates of p ' and p are 

ty+i© 


(2a) 

(2b) 


p ' = e 


Wn+i0 

p = e 

Equation (2a) transforms the semi- infinite periodic strip in the 
z-plane, bounded by 0 = 0, 9 = 2ji, and the corrugation ty(0), 

into the entire region outside the p' -plane near circle that corre- 
sponds to the corrugation. The corrugation given by the Cartesian 
coordinates and 9 in the z-plane is represented in the p‘ -plane 
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by the near circle with polar coordinates and 0. Similarly, 

equation (2b) transforms the semi- infinite periodic strip, bounded 
by 0 = 0, 0 « 2tc, and the 0-axis (^q = 0), into the entire 

region outside a unit circle in the p-plane. The CMF z - £ becomes 
a function that is regular everywhere outside the p-plane unit 
circle and is therefore expressible by an inverse power series . 


z 


t = log (p'/p) 



(3) 


wherein the constant term, representing a relative translation 
between the z-plane and the £ -plane, has been made zero. The trans- 
formation log (p'/p) in equation (3) has been used in reference 2. 
On the boundaries equation (3) becomes, with p = e 1 ^ and 



OO 


CO 


W) 

= y / a n 

cos n0 

+ 

sin n 0 


1 


1 



OO 


, CO 



•£> n 

cos n0 

- > ' a n 
<C r n 

sin n0 


1 1 


(4) 


The mapping fun.ction '4/(0) - i£(0) for a given boundary 'fy(Q) 
can be obtained from equations (4) . Conversely, special families 
of corrugations are obtained by selecting various harmonics in equa- 
tions (4); for example, a simple type of corrugation is given by 


\J/(0) = _ £•[ cos 0 
2 

£(0) = - sin 0 > 

9(0) = 0 - e(0) = 0 + ~ sin 0 


(5) 


where T/2 i 3 the thickness ratio of the corrugation, defined as 
the total height h (fig, 1) divided by the wave length 2it| the 
quantity T is thus analogous to the thickness ratio 
(maximum thickness/chord) of airfoil sections. The members of this 
family (equations (5)) corresponding to T=0, 0.1, 0.2, and 0.3 


are plotted in figure 2 as 


* 

jcT/2 


against 0 . 
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Once the mapping function of a corrugation is known, the veloc- 
ity v at any point on the surface is given by the product of the 
velocity on the straight-line boundary, which is unity, and the 
stretching factor jd£/dz| from the straight line to the corruga- 
tion; thus, by vise of equations (l). 


v 


dXl 1 

dz 1 1 dil/ + id0 1 

I lift | 



For the special family of corrugations given by 
the velocity distribution reduces to 


( 6 ) 


equations (5) 


Av _ v-1 
jxT/ 2 “ jtT/2 


2_ 

r i 

JtT 

i . . ' : 

/ 

A 1 



1/1 t 

3T ) + rtT cos 0 


H 

\2 } 


1 


( 7 ) 


Figure 3 shows the velocity distributions of members of the special 
family shown in figure 2 . As 1 — > 0 , and Av also — > 0 , but 
Av 


■0/ 

am 


and 


X.TJ2 


4— cos 9. These limiting values agree with 


the thin-airfoil results obtained by Allen (reference 1). 


In the general case of a given arbitrary corrugation ( 9 ) , 

the GMF il/(7i) — ie(0) can be determined 1 'oy successive approxima- 
tions. Suppose, for example, that the zeroth approximation to the 
corrugation ^'(S) is the straight line 'l' o (0) — i£ Q (0) = 0. The 

first-approximation ordinates ^(0), corresponding to a set of 
evenly spaced 0 values, are then obtained from the given boundary 
at the abscissas 0 q = 0. The function £^(0), conjugate to 
* 4^1 ( 0 ) , is determined by harmonic analysis and synthesis in accord- 
ance with equations (4), The resulting first-approximation 
CMF — i€j_(0) yields the coordinates \j/-j_(0), 0^(0) =0—^(0) 

of a boundary, which is compared with the given boundary. If the 
agreement is not satisfactorily close, the procedure is repeated; 
the second-approximation ordinates ^(0) corresponding to the same 
set of evenly spaced 0 values, are obtained from the given boundary 
at the abscissas 0^(0) ~ 0 — e i(0)^ etc. 
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As an example of the general procedure and for comparison with 
the results of Allen (reference 1) in the case of nonnegligible thick- 
ness , a cosine corrugation with a thickness ratio of 0.2 was taken 
as )jr = — ~ cos 6. The zeroth approximation was chosen as the 
special corrugation T = 0.2 cf figure 2. The maximum difference 
in ordinates of the two surfaces was 30 percent. After the first 
approximation the maximum difference between the ordinates of the 
given cosine corrugation and the first -approximation boundary was 
reduced to about 4 percent. A second and a third approximation 
further reduced the difference to 1 and 0.25 percents , respectively. 
The resulting CMF and the velocity distribution for the third 
approximation are given in table 1. Figure 4 shows this velocity 
distribution as well as the approximate velocity distribution based 
on thin-airfoil theory. A.s was demonstrated by Allen, the approxi- 
mate velocity distribution is a cosine distribution. 

The results obtained by the two methods for a thickness ratio 
of 20 percent differ appreciably; the maximum difference is about 
16 percent of the maximum increment of velocity over the free- stream 
value. In the range of thickness ratios contemplated by Allen, how- 
ever, the results of thin-airfoil theory are undoubtedly of sufficient 
accuracy, as far as incompressible potential flow is concerned. 


TEED FLOW 0WR A BUMP 


Consider a surface that is perfectly flat except for an isolated 
bump or a disturbance of constant chord length and infinite span; 
assume the flow over the surface to be at right angles to the span 
and of magnitude unity sufficiently far from the bump. A two- 
dimensional symmetrical flow is obtained by reflecting the bump in 
the plane surface. This problem is solved by conformally mapping 
the symmetrical section, taken in the z-plane, into its axis of 
symmetry, taken in the £ -plane. The coordinates in the two planes 
are (fig. 5) 


z = x + iy 

i = % + ^ 


( 8 ) 


The CMF z - £ becomes a function regular in the exterior of a 
circle jp| = R as a result of the Joukowski transformation 


t 
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and can "be expressed as an inverse power series 

as 

z-£ = 4X + iiy=2^ 


0 


(9) 


For corresponding points on the "boundaries in the z~ and 

p -plane s. 


ix w = 2 ^ oos ^ + 2 ^ sin ^ 


o 


CO 




• = 2 s oos ^ - 2 55 sln nl t 


E n 


E n 


( 10 ) 


x(0) = r cos 0 + Ax(0) (r = 2R) 


The velocity v at any point of the symmetrical section is the 
product of the velocity at the corresponding point of the circle 
jp| = E and the stretching factor jdp/asj . The result is 



The sections under consideration are now assumed to he symmetrical 
with respect to "both the coordinate axes and to have a horizontal 
tangent at their chordwise extremities on the x-axis. The Fourier 
series (equations (10)) are thereby simplified; symmetry with 
respect to the x-axis requires that b n = 0, and symmetry with 
respect to the y-axis requires the vanishing of even harmonics. 
Hence , 




Go 

vn Q rx a 

Ay = - > — sin n0 

H n r 

1 


> (n odd) 


(12a) 

(121>) 
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The condition that the section have a horizontal tangent at the 
chordwlse extremities, that is, 

V naQ cos n0 

dy _ 1 

dx a? 

r sin j) + V na n sin n0 

1 

he zero for 0 = 0, is satisfied if 


CO 00 

2 na n = 0, r +2 n 2 a n 0 0 (13) 

1 1 

A simple example of a family of humps satisfying conditions (13) 
is given ty 



x = (l - cos 0 + | cos 30 (15) 


where the value of r has been so adjusted that the chordwise 
extremities of the section are at x = +1. The thickness ratio T 
is defined as twice the height of the hump divided hy its length 2, 

The humps given hy equations (14) and (15) are shown in figure 6(a) 
for T = 0, 0.1, 0.2, and 0.3. The family of symmetrical sections 
shown in this figure was derived hy Kaplan in reference 5 hy a 
generalization of the Joukowski transformation. The corresponding 
velocity distributions are shown in figure 6(h). 

The velocity distribution on the wall (that is, for |x| >1, 
y = 0) can he obtained from the general expression for the velocity 
at any point in the plane outside the section,* thus, for an arbi- 
trary airfoil situated at an angle of attack a in a free stream 
of unit velocity, the expression for the derivative w z of the poten- 
tial function in the airfoil z-plane is 

Wz = dz/dp 



I 


where Wp is the derivative of the potential function in the circle p-plane. If p . is 
witten in the form 

p = 

i 


0 = V - Hr 0 J 

evaluation of Vp and dz/ dp yields the following formulas for the magnitude vz and the 
direction n of the velocity vector in the airfoil plane: 



where {3^ is the zero-lift angle. 

For a symmetrical flow and section, a = (3j = 0; whereas, for points on the wall, 
0 = ^Az/^0 = 0. Equation (17) thus reduces to 
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(v z ) 


wall 


sinh Q 
sirAi Q + 


(19) 


The value of f^~r) 


for arbitrary Q>>0 is obtained by differen- 


ce) 


tiating equation (12b), replacing B with e^, and using for a n 
the values .previously determined for the section. 


The conformal transformation of a sinusoidal bump was next 
determined by the method of successive approximations outlined in 
the preceding section. The thickness ratio T was taken to be 
0.2, so that the symmetrical section has the equation 

y = ±0.1 (1 + cos it x) (-1< X < 1) (20) 

The symmetrical section of figure 6(a) with T = 0.2 was chosen as 
the zeroth approximation. Two approximations were carried out. The 
maximum differences between the ordinates of the given boundary 
(equation (20)) and the successively derived boundaries were about 
5, 1.5, and 0.25 percent for the zeroth, the first, and the second 
approximations, respectively. Table 2 contains the data for the 
second approximation and figure 7 shows the velocity distribution 
over bump and wall. The approximate result of Allen, obtained on 
the basis of thin-airfoil theory, is also shown in figure 7. The 
maximum difference between the two curves is about 8 percent of the 
maximum increment of velocity over the free- stream value. This 
difference, it should be remembered, is for a 20-percent thic]giess 
ratio; for the very small thickness ratios considered by Allen, 
thin-airfoil theory is quite adequate. 


Figures 6 and 7 show that the velocity distribution in the 
neighborhood of the point where bump meets wall merits discussion. 

A symmetrical section can become tangent to the wall in any one of 
three ways: with infinite curvature, zero curvature, and finite 

nonzero curvature. The case of infinite curvature, properly called 
a cusp, holds for the special family of sections given by equa- 
tions (14) and (15) and also for the trailing edge of a symmetrical 
Joukowski airfoil. The velocity curve corresponding to both bump 
and wall has a minimum value at the cusp; the velocity gradient at 
the cusp is finite on the bump side and infinite on the wall side 
as indicated in figure 6(b). See reference 4 for a comparison with 
experiment . 


Zero curvature is obtained at the sharp edge of the symmetrical 
section if, in addition to conditions (13), the following equation 
holds ; 
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The velocity curve in this case has a minimum at seme poinb on the 
hump and the velocity gradient is continuous at the point where the 
section meet3 the wall. 

The section obtained by reflecting the cosine bump in the wail 
has finite nonzero curvature at its sharp edge. The velocity is a 
minimum at a point on the bump, as seen in figure 7 , but the conti- 
nuity of the velocity gradient at the sharp edge is, from the cal- 
culations of this paper, still an open question. It is conjectured 
that , at a sharp edge of this type, the velocity gradient has a 
finite discontinuity. 


EXTERIOR- IETSEI0R BUMPS 

The flow over a bump has been derived in the preceding section 
by reflection of the bump contour in the wall and analysis of the 
resulting symmetrical section. If an ex berior- interior bump, namely, 
a distortion of part of a wall in both directions perpendicular to 
the wall , is reflected in the wall, a symmetrical figure- eight sec- 
tion results, as indicated in figure S. The mapping of such a con- 
tour onto a circle can be accomplished as previously described. 

It appears, hoi-fever, that the derivative dz/dp of the transforma- 
tion will be zero at a point outside the circle corresponding to a 
point within the loop consisting of the interior part of the b’mip 
contour and its reflection. That such a zero must exist becomes 
evident upon tracing the paths around the figure-eight contour 
corresponding to concentric circles larger than the basic circle. 

As Indicated schematically in figure 3, the transition contour 
between those of figure- eighu type and those simply connected has a 
sharp- edged extremity at the point F inside the loop formed by the 
interior part of the bump contour and its reflection. At this 
sharp edge, dz/dp = C. Although this property of a looped contour 
might be useful, for example, in locating the singularities of a 
mapping function, this method of attach does not yield the desired 
flow over the exterior- interior bump (the flow actually obtained is 
that whose zero streamline Is the path ABCDEFGE In fig. S). It may 
be noted that the conventional application of thin-airfoil theory 
also breaks down in this case. 


The flow over an exterior- interior bump can be obtained by 
mapping the bump and the wall contours onto an infinitely long 
straight line. Points on the bump contour are related to points on 
the straight line by the CMF, 
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z - £ = Ax + iA y 

a3 indicated in figure 9, and points on the straight line are related 
to points on the unit circle by the bilinear transformation, which 
for the problem under consideration is taken as 



( 21 ) 


Equation (21) transforms the upper half £ -plane into the region 
exterior to the unit circle ; corresponding points of both regions 
are shovm in figure S. For points on the boundaries, the inyerse- 
povrer- series expression for z - £ (equations (9) and (21)) yields 


Ax 

CP 

= 2 

a n 

cos n0 

CO 

sin n0 


0 



0 



CD 



co 


Ay 

_ *\' i 
~ 2L i 

^n 

cos n0 

-2*n 

sin n0 


0 



0 


X 

= I + Ax 

= tan 

90° - 0 
2 

+ Ax(0) 


y = Ay 






( 22 ) 


The Telocity distribution at the surface of the bump is obtained 
from the complex Telocity function wz: 


\r 


v z = 


L. - 


dz 


it ^ 1 - & • •- & * » 


4AA 
aj<i ) 


The absolute magnitude of the Telocity t is therefore glren by 


T = 


1 - (1 + sin 0) 


dAx 

djZi J 


l2 


(1 + sin 0) , 

j 


As a simple example of an exterior- inter ior bump, the family 
represented by 
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Ax = b sin f> 

Ay = y = b cos 0 

90° - 0 v , hc 
x = tan g — ^ + b sin p 

Is illustrated for various values of b, together -with the corre- 
sponding velocity distributions, in figures 10 and 11. By the 
neohods previously described and also by superposition of solutions 
by linear combinations of CME's, arbitrary distortions of a straight 
nail may be analyzed or synthesized. 


COI'TCLIBIOh 


The velocity distributions on corrugations and bumps as deter- 
mined by conformal transformation are, in the case of 20-percent 
thickness ratio, appreciably different from the corresponding results 
by thin-airf o'j 1 theory. The maximum differences, expressed a3 frac- 
tions of the maximum increment of velocity over free- stream velocity 
produced by the disturbance, amount to approximately IS percent for 
a sinusoidal corrugation of 2C-percent thickness ratio and 3 percent 
for a, sinusoidal bump of the same thickness ratio. In the limit of 
zero thickness ratio, the results by conformal transformation are 
identical trlth the results by thin-airf oil theory. 


Aircraft Engine Besearch laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 


EEESREPTCSS 

1. Allen, H. Julian: Botes on the Effect of Surface Distortions on 

the Drag and Critical Mach E umber of Airfoils. 1JACA ACE 
Bo. 3129, 1943. 

2. Theodorsen, Theodore: Theory of Wing Sections of Arbitrary Shape. 

iJACA Bep. Bo. 411, 1951. 



MCA ABE Ho. E5A05 


13 


3. Kaplan, C.f The Flow of a Compress ill e Fluid Past a Curved Sur 

face. BACA ABE Bo. 3K02 , 1943. 

4. Preston, J. 3., and Sweeting, H. E.: The Experimental Determi- 

nation of the Boundary Layer and Bake Characteristics of a 
Simple Joukowski Aerofoil, with Particular Deference to the 
Trailing Edge Legion. F. M. 570, S. & C. 1499, Ae. 2168, 
(British), March 1943. 



14 


NACA AER No. E5A05 


TA3IE 1, - CMF AND VELOCITY DISTRIBUTION FOR COSINE CORRUGATION 
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TABLE 2. - CMF AND VELOCITY DISTRIBUTION FOR COSBIE BUMP 
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Figure ( - Conformal transformation of a corrugated surface . 
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Figure 2.- Special family of corrugations, 
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Fig. 3 



Figure 3.- Velocity distributions of special 
family of corrugations. 
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Fig. 5 



Figure 5. — Conformal transformation of a symmetrica! section. 
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Fig® 6 b 







NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 



NACA ARR NO. E5A05 










NACA ARR No. E5A05 


Fig. 9 



Figure 9.— Mapping of bump using b/tinear 

transformation . 
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Figure ll.« Velocity distributions of special family of infinitely extended exterior 
interior bumps « 
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